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Abstract: Scanning specimens in liquids using commercial atomic force microscopy 
(AFM) is very time-consuming due to the necessary try-and-error iteration for determining 
appropriate triggering frequencies and probes. In addition, the iteration easily contaminates 
the AFM tip and damages the samples, which consumes probes. One reason for this could be 
inaccuracy in the resonant frequency in the feedback system setup. This paper proposes a 
frequency function which varies with the tip-sample separation, and it helps to improve the 
frequency shift in the current feedback system of commercial AFMs. The frequency 
function is a closed-form equation, which allows for easy calculation, as confirmed by 
experimental data. It comprises three physical effects: the quasi-static equilibrium condition, 
the atomic forces gradient effect, and hydrodynamic load effect. While each of these has 
previously been developed in separate studies, this is the first time their combination has 
been used to represent the complete frequency phenomenon. To avoid "jump to contact" 
issues, experiments often use probes with relatively stiffer cantilevers, which inevitably 
reduce the force sensitivity in sensing low atomic forces. The proposed frequency function 
can also predict jump to contact behavior and, thus, the probe sensitivity could be increased 
and soft probes could be widely used. Additionally, various tip height behaviors coupling 
with the atomic forces gradient and hydrodynamic effects are discussed in the context of 
carbon nanotube probes. 

Keywords: frequency shift function; jump to contact; liquid environment; atomic 
force microscopy 



Sensors 2014, 14 



9370 



1. Introduction 

In the last ten years, the atomic-scale resolution of frequency-modulation atomic force microscopy 
(FM AFM) has been taken advantage of to image and conduct force spectroscopy measurements of 
biological samples in liquids. This has enabled quantitative measurements of conservative and 
dissipative forces, and subnanometer resolution imaging with piconewton-order loading forces at the 
solid and liquid interface. The thermal noise in the liquid environment increases the vertical vibration 
amplitude of the probe, hindering the FM AFM techniques from adopting probes with stiffness greater 
than 40 N/m [1]. The stiff probe is sensitive enough for strong interactions; however, weak interactions 
require a soft probe. Soft probes are highly susceptible to the problem of "jump to contact," which is 
minor reason why the soft probe has not been adopted in FM AFM. On the other hand, Sader et al. [2] 
successfully overcame the amplitude problem of soft probes and constrained the vibration amplitude due 
to thermal noise under 2 nm for a probe stiffness of 3 N/m. This achievement indicates that the thermal 
noise problem of soft probes is solvable. To overcome the instability of "jump to contact", this paper 
provides a function which can calculate frequency shifts, thus predicting "jump to contact" occurrences. 
The proposed frequency function could also be used in the tapping-modulation (TM) AFM and other 
forms of AFM. This frequency function is important for feedback control of AFM, and it can be used for 
the inverse calculation of tip-sample interaction forces in liquids. 

FM AFM is a promising technique to measure biological samples. Recently, Bruker [3] developed a 
new TM AFM technique, referred to as Peak Force Quantitative Nanomechanical Mapping (PF QNM), 
applied to measure living and soft samples. PF QNM, succeeding from TM AFM, replaces detecting 
frequency (or phases) shifts with the peak force as the measurement feedback. Different from typical 
AFMs driven around resonance frequencies (around a hundred kilohertz), PF QNM is driven at a 
frequency far below (a few kilohertz). The probe amplitude of PF QNM is around hundred nanometers 
like general TM AFM, but that of FM AFM is around few nanometers. PF QNM features in measuring 
wide range of Young's moduli of biological samples with noise lower than general TM AFM. However, 
FM AFM, because of its small tip amplitudes of tip, has subnanometer resolution for measuring 
biological samples [1]. 

A probe selected in biological AFM depends on modulations of AFM and the resolution required for 
the specimen. Most atomic force microscopy surface imaging of bio-specimens is performed by the 
contact-modulated approach. To avoid destruction of the biological specimen [4], probe stiffness 
remains within the low range, 0.03-0.5 N/m [5,6]. However, low-stiffness probes easily jump to contact 
and suffer impaired resolution. For this reason, interest in the use of noncontact AFM has grown. FM 
AFM, one type of noncontact technique, uses low amplitude and measures the frequency shift of the 
probe, enabling higher resolution than other techniques. Within a liquid environment, atomic resolution 
can be achieved with biological applications of FM AFM, achieving a maximum signal by initial values 
of amplitudes at 0.6-1.0 nm and leading a resonance frequency shift of tens of hertz [1]. Thermal noise 
makes control of the amplitude in the liquid difficult; the peak value of the power spectrum density of the 
thermal vibration is given by (z t h) P eak = v /2k B T Q/ Tr fok, where fo and Q are resonance frequency and the Q 
factor of the microcantilever, k B and Tare Boltzmann's constant and absolute temperature, and k is the 
probe stiffness. Typically, we increase stiffness (the resonance frequency) to decrease the thermal 
amplitude of probe, thus, avoiding "jump to contact" issues. However, increasing stiffness decreases the 
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sensitivity when measuring the atomic force. Reversing the assumption on increasing stiffness, we try to 
reduce the value of Q (for example, by improving the microcantilever shape) and consider the benefit of 
"jump to contact." It may be possible to develop a highly sensitive probe in liquid. Thus, this 
paper adapts the use of soft probes and provides a frequency shift function which predicts "jump to 
contact" behavior. 

The frequency function setup in commercial AFM devices can be improved easily. Scanning 
specimens in liquids using commercial atomic force microscopy (AFM) is time-consuming due the need 
for a try-and-error iteration to determine appropriate triggering frequencies. In addition, this iteration 
easily contaminates the AFM tip and damages samples, which consumes probes. One significant reason 
for this is that the frequency function is not accurate. The frequency function is not a function of the 
separation between tip and specimen. This results in the TM AFM incorrectly using resonant frequency 
to trigger the microcantilever and recording the wrong signals to judge height of specimen when the tip 
approaches the sample. The resonance frequency of AFM probes in liquid has been commonly 

2 1/2 

approximated by f=f mc (1 + npjb I Am) , where f vac is the resonance frequency in vacuum, pf and m are 
the fluid density and mass per unit length, and b is the width of the probe. This approximation is suitable 
for the probe in far field, but it does not adequately represent the resonance frequency when the probe 
approaches close to the specimen. Unfortunately, this equation is used in commercial AFM devices, and 
the differences could be on the order of a few hundred hertz. Two causes underlie the shift in resonance 
frequency in liquids: (1) the interaction of atomic forces between the tip and the specimen, and (2) the 
disturbance of liquid in the space between the probe and the boundary surface on which the specimen is 
placed. In the first, atomic forces manifest over a few nanometers, and the effects to the shift can again 
be divided into two types: (a) the probe's eigenfunction varying with the tip-sample separation (the 
quasi-static equilibrium condition [7]), and (b) the force gradient effect around the amplitude region 
(Giessibl's theorem [8]). The effect of the disturbance of liquid is mainly based on Green and Sader's 
theorem [9] and arises from the liquid pressure on the probe varying with probe-surface separation. This 
effect manifests over tens of micrometers. Sader et al. suggested a frequency equation that considers the 
atomic force effect and the liquid pressure effect [1]. The aforementioned frequency function contains an 
opened-integration (j z ), which is ultimately problematic for AFM engineers. To improve it and allow 

for the use of soft probes, this paper provides an equation for the frequency shift as a function of 
tip-sample separation. The frequency function is a closed-from equation and is quick to obtain results. 
Experimental data are also proven with the proposed frequency function. Additionally, various tip height 
behaviors coupling with the atomic forces gradient and hydrodynamic effects are discussed in the 
context of carbon nanotube probes. 

2. Mathematic Model and Discussion 

2.1. Quasi-Static Equilibrium Condition 

The equation of motion from the Euler-Bernoulli beam is: 



dx 



j[EI(x) 



d 2 w(x,t) 
dx 2 



] + m(x) 



d 2 w(x,t) 
~~di 2 ~~~ 



+ c(x) 



dw(x, t) 
dt 



= p(x,t) 



(1) 
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The following discussion is limited to beams with uniform properties along their lengths, i.e., 
Young's modulus, moment of inertia, mass per unit length, and the damping coefficient given by the 
constants: E, I, m, and c. The total displacement is w(x,t) as the sum of displacements that would be 
induced by static application of the support motion, i. e. , the quasi-static displacement, plus the additional 
displacement due to the dynamic inertial and viscous force effects [10]. At the clamped-end, 
w(0,t) = 0 and w'(0,t) = 0. At the tip-end, affected by the atomic force, w"(0,f) = 0 and 
EIw"'(L,t) = m tip <k(L,t) + k ts (z)w(L,t), in which m tip is the mass of the tip, k ts (z) is the gradient of the 
force curve, and L is beam length. Applying the above boundary conditions to Equation (1) gives the 
eigen function: 

(3R k - R ti X )(sin p n cosh fi n - cos (5 n sinh (5 n ) + (1 + cos /?„ cosh fi n ) = 0 (2) 

where /?„ is a coefficient in the separation of variable method, and the value is relative to the natural 
frequency of the microcantilever. The stiffness ratio and the tip-mass ratio are Rk = —k ts /ko and 
Rtip = nttiplmL. ko = 3EI/L 3 is the stiffness of the microcantilever. If k ts {z) is assumed to be specific form 
(for example, Lennard- Jones form), fi n can be solved by numerical method. The eigenfrequency is 
defined by f n (R k ,R up ) = J^(p n /L) 2 /2-rt. The first three modes of the eigenvalues fii, and versus the 

stiffness ratio for certain tip-mass ratios are plotted in Figure 1 . All eigenvalues are normalized by /?o, the 
eigenvalue when Rk = 0. Thus fto leads to fo, which represents the resonance frequency in the far field. 
From above, the resonance frequency is not constant, as it varies with tip-sample separation. 
Experimental data obtained from Gotsmann et al. is also plotted [11]. The minimum stiffness ratios in 
the associated experiments are shown to be a little bit less than zero (Rk — > 0 ). 

Figure 1. The first three normalized eigenvalues /?„//?o varying with stiffness ratio Rk and 
tip-mass ratio Ru p . Rk > 0 represents a contact, -1 < Rk < 0 represents a oscillation, and 
Rk<—\ represents a jump. Experimental data, denoted by *, are from Gotsmann et al. [11]. 




Stiffness ratio R, 



However, the first eigenvalue may sometimes have a trivial solution wheni^ < — 1. A critical point, 
(Rk, fix) = (~1,0) is utilized to determine where static equilibrium is broken and a jump occurs. This can 
be proven by applying an asymptotic series in Here, substituting sin/?i = /?i - /?i 3 /3! + 
cos/?i =fii- y?i 2 /2! + -, sin/z^i = fix + ^i 3 /3! + and cos/z/?i = + fi\l2\ + - into Equation (2) gives 
(3Rk- Rti P P\ 4 )(2Pil3) + 2fi\ = 0, in which the higher order terms are neglected. Setting /?i — > 0 leads to 
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Rk = — 1. This implies that the absolute value of the gradient of atomic force is the same as the stiffness 
of the beam. This corresponds to the case where the tip approaches the sample from afar, and the 
stiffness ratio decreases from 0 to -1. The minimum stiffness ratio of a stiff microcantilever may not 
reach -1, but a soft microcantilever may reach -1. Then a small disturbance brings the soft cantilever 
into the region Rk< -1, where the attractive force exceeds the restoring force, accelerating the tip toward 
the sample. In Figure 1, this nonlinear result breaks the assumption that the frequency shift is a linear 
function of the tip-sample gradient in the range [— 1, 0], Af = f 0 /2(k r Jko) as provided [12]. Figure 1 also 
shows that the second and third eigenvalues remain constant with the stiffness ratio. It is also clear that 
the tip-mass ratio influences the eigenvalues, but is independent of the tip-sample separation. From an 
experimental data point of view, typical probes are conservatively designed according to (Rk — > 0~). 
Since the performance of the probe can be predicted in our proposal, the stiffness ratio range for an 
effective probe can be expanded to -1 < Rk< 0. 

2.2. Atomic Force Gradient Effect 

Moreover, the frequency is a function of vibration amplitude, especially oscillation in the near field 
featuring rapid variation of the atomic force. Giessibl modified the results of Albrecht et al. [12] and 
utilized canonical perturbation theory to solve the frequency shift for large amplitudes [8]. His solution 
assumes that the gradient of the restoring force is larger than that of the atomic force, and that the 
cantilever motion is well described as approximately harmonic. Then the frequency shift is strictly 
proportional to^ and l/ko, i.e., Af = (f n /2ko)<k*(A 0 )>, in which the bracket, < >, indicates averaging 
across one oscillation cycle and is the weighted average of the gradient of the force: 

(** W) = j J o K (x)^A 2 0 -x 2 dx (3) 

Giessibl applied to the atomic resolution of his studies, but the equilibrium of the quasi-static 
conditions was not considered [13—15]. To remains constant in his papers. However, his experimental 
frequency shifts were still well matched with the theoretical prediction. Because the probes employed in 
these experiments were sufficiently stiff, Rk — > 0~, so f was constant. In this paper, the combination of 
eigen frequencies obtained from the quasi-static equilibrium and Giessibl's theory leads to the resonance 
frequency in vacuum: 

/«(^V^)=/-(^VP + -^Jl^^-W4-^] (4) 

If k rs (z) is assumed to be Lennard- Jones' formula, the integration could be solved to be a function of 
amplitude. Figure 2a shows the frequency shift of a tungsten tip on a KC1 (100) surface in which 
Equation (4) is applied to compare with experimental data from Giessibl's study [16]. Model data was 
generated with an amplitude A 0 = 0.15 nm, a stiffness of ko = 1800 N/m, length 2.4 mm, width 130 |un, 
thickness 214 jam, and eigenfrequency set to f = 25,068 Hz. 
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2.3. Hydrodynamic Effect 

The modeled frequency shift takes into account hydrodynamic effects when the probe is immersed in 
liquid. Sader [9] extended the boundary integral technique of Tuck [17] and presented an explicit 
semi-analytical theory to solve for an oscillating cantilever immersed in a viscous fluid nearby the 
specimen surface. This paper extended Sader's work to FM AFM probe. It assumes that the fluid was 
incompressible and the oscillation amplitude was small. The fluid flow around the microcantilever is 
governed by the incompressible unsteady Stokes equation: 

p f ^- = -VP + r J V 2 v (5) 

where v, P, and rj are vector velocity, pressure, and viscosity of the fluid, respectively; and the Reynolds 
number of the flow is R e = npfhb I2n [18], and the Reynolds number of a soft AFM probe of width 
30 jam oscillating in water is in range [0.1, 20]. Fluid velocities are restricted to be normal to the solid 
surface, whereas no-slip boundary conditions are enforced at interfaces. The hydrodynamic load 
per unit is obtained from pressure differences between the top and bottom of the beam, 
Ph = nw(x,t) J^ 2 /2 AP(Od^, in which w(x,t) is the normal velocity and AP(<f) is the pressure jump 
across the beam. Let the Fourier transform be X = f°° Xe lwt dt. Taking the Fourier transform of 

J — oo c 

Equation (1) yields EIw - marw + icocw = p. p is the hydrodynamic loading component. Examining 
the Fourier transformed equations of Equation (5) and setting V • v = 0 give the general form of the 
hydrodynamic loading component as p = n 2 f 2 p f b 2 T(Re,b,h)w, where h is the gap between the probe and 

the surface and T(R e ,b,h) = ij^ 2 /2 AP(Qd^/(nR e ) is the complex hydrodynamic function [19]. Its real 

part T r represents the effect of added mass on the surrounding fluid, and its imaginary part T, represents the 
damping effect of the fluid. The hydrodynamic load equation is formally exact in the limit of 
0 < Lib « 1. Substituting p into the Fourier-transformed motion equation leads 

2 — 1/2 

fh(R e ,b,h) = fac [1 + npjb Tr(R e ,b,h)/4m] . Note the other formula for the frequency function can refer 
to the study by Naik et al. [20]. Thus f vac replaced by Equation (4) gives the frequency: 

1 + 1 k ts (z - x)^A 2 -x 2 dx 
f h (R k ,R tip ,A 0 ,R Q ,b,h) = f n (R k , R tip ) * o4 ° b2 (6) 

[1 + ^— r r (Re,M)] 1/2 
4m 

This is the desired equation, representing frequency as a function of the eigenfrequency and 
proportional to the atomic effect and inversely proportional to the hydrodynamic effect. From an 
engineering point of view, the interaction force curve can be assumed and the probe properties and the 
amplitude of the probe are given by users. This frequency function works while the hydrodynamic 
function Y r is given. Here, Y r has an approximation function and will be discussed later. Note that T r is a 
function of the Reynolds number R e , which itself is a function off,. Applying the numerical method 
allows to be easily solved. As a result, Equation (6) is useful in commercial devices. This equation is 
different to the equation provided by Sader et al. [1]. It contains the quasi-static equilibrium condition 
and can predict the performance of the soft AFM probes. Furthermore, the inverse method for 
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calculating the tip-sample forces in liquid from the experimental frequency shifts can be realized by 
applying Giessibl's technique [16] along with Equation (6). 

2.4. Hydrodynamic Function 

The details of the hydrodynamic function T r (R e ,b,h) can be obtained from the works of Green, Sader, 
and van Eysden [9,19] who provided closed-form analytical expressions for the entries of all 
submatrices except four complicated terms, {A\^2,C\,C^). In this paper, these four terms were computed 
by applying numerical Gauss-Legendre quadrature. But it was a time consuming procedure. 
Tung et al. [21] calculated the hydrodynamic function and provided the polynomial approximation 
function by curve fitting in the domain R e e [ 1 0~ 2 , 1 0 4 ] . These fitting results were found to have errors on 
the order of a few percent, which might shift the frequency a few tens of hertz (the eigenfrequency is 
around one hundred kilohertz). 

Figure 2. Frequency shift versus separation for: (a) a vacuum case with prediction by the 
proposed method (solid line) and experimental data (dots) [16], and (b) a liquid case with 
prediction by proposed method (solid line), experimental data (dots) [2], and results (dashed 
line) by the proposed method but with Y r replaced by approximation [21]. 
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Figure 2b shows the frequency shift of a microcantilever immersed in water plotted against the gap 
between the probe and the surface. The dimensions of the microcantilever (Nanosensors EFM 
cantilever) are length 225 \im, width 23 [im, thickness 3 jam, a tip height of 13 jam, a oscillation 
amplitude of 1 nm, and the resonance frequency in water is 13.14 kHz. The solid line plots results 
obtained from Equation (6) directly by the semi-analytical approach [9]. The dashed line represents 
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results from Equation (6) but with the hydrodynamic function derived from Tung's polynomial 
approximation [21]. The dots are experimental data obtained from Sader et al. [2]. The dashed line 
deviates when the gap is larger than 10 (J,m because Tung's approximation is outside the convergent 
region. Tung's approximation saves time and is useful in predicting the shift tendency of the resonance 
frequency. When the gap is much smaller than the microcantilever width, a change in the gap has less of 
an effect on the resonance frequency due to the hydrodynamic effect. Comparing Figure 2a with 2b, the 
atomic force affects the frequency shift by several hertz in a region a few nanometers in the near-field, 
and the hydrodynamic affects the frequency shift by hundreds of hertz in the region around a hundred 
micrometers. As a result, while the tip height is a few micrometers, these two effects are coupled. 

2.5. Tip Heights and Frequency Shifts 

Figure 3a shows frequency shift versus tip-sample separation for various tip heights, set at 0.1, 1, and 
10 (j,m. A case of a silicon tip and a polystyrene surface is studied with probe dimensions of length 
250 jam, width 35 jam, and thickness 5.7 um. 

Figure 3. Frequency versus tip-sample separation for (a) probes with various tip heights, 
and (b) a CNT-tip with various microcantilever lengths. 
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The resonance frequency in water is 232.5 kHz. The Lennard- Jones interaction force is assumed to be 
F(z) = Jji?/180z 8 -A 2 R/6z 2 , in which R = 20 nm, A\ and A 2 are the Hamaker constants, 
0.838873 x 10" 70 J m 6 and 1.15072 x 10" 19 J, respectively [22]. The results indicate that the h tip = 0.1 \xm 
case presents the hydrodynamic affects from far field to the near field, tip-sample separation between 
0.1-10 Lim. Hydrodynamic effects are effectively constant when the tip-sample separation is less than 
0.1 (j,m. Moreover, typical tips are around ten micrometers, and the frequency shifts due to the 
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hydrodynamic force are on the order of hundreds of hertz, as shown h t i p =10 jam. Hence, high tips 
feature reducing of the hydrodynamic and presenting the atomic force effects. 

Recently, tips extended by a carbon nanotube (CNT) positioned on the tip have been used to measure 
high- aspect ratio samples. The CNT can be hundreds of nanometers to a few micrometers long, and a 
few hundred nanometers long is suggested to avoid buckling. In the following calculation, the total 
height of a CNT-tip is set to 10 jam, and the rest of the probe geometry is the same as in previous 
calculations. The carbon-carbon interaction potential is F(z) - 4 e^ncr 2 [-(a / z) 10 /5 + (o7z) 4 / 2] , 

where the parameters are assumed to be s = 4.751 xlO" 22 eV, n = 0.114 A" 3 , and o = 3.407 A [23]. To 
study the jumps, the microcantilever length was set sequentially to L = 125, 250, and 500 [im (ko 
distribution ratio is 1/8:1:8) as shown in Figure 3b. Since the CNT-tip has a rapid interaction, the 
stiffness ratio becomes less than -1 in the first two cases, and that represents CNT-tip jumps. The 
frequency shifts and critical points of the jumps can also be analytically obtained by the proposed method. 

3. Conclusions 

A frequency function for AFM in liquid has been introduced, with the frequency shift derived as a 
function of six parameters: stiffness ratio, mass ratio, oscillation amplitude, Reynolds number, probe 
width, and tip height. This proposed frequency function is in closed form, which is easy for engineers to 
implement in AFM devices. It can predict the jump to contact behavior. Accordingly, probes with 
stiffness ratios smaller than -1 could be applied; that is soft probes could be easily used to increase the 
probe sensitivity. By application of the frequency function, the probe with high tips demonstrates the 
effect of coupling between atomic and hydrodynamic forces, and soft probes with CNT-tips show 
jumping behavior in accord with the predictions. 

Acknowledgment 

The author wishes to thank National Science Council of the Republic of China, Taiwan, for financial 
support under contract NSC- 1 0 1 -222 1 -E-390-003 . 

Conflicts of Interest 

The authors declare no conflict of interest. 
Reference 

1 . Fukuma, T.; Jarvis, S.P. Biological applications of FM-AFM in liquid environment. In Noncontact 
Atomic Force Microscopy; Morita, S., Giessibl, F.J., Wiesendanger, R., Eds.; Springer- Verlag: 
Berlini/Heidelberg, Germany, 2009; pp. 329-345. 

2. Sader, J.E.; Uchihashi, T.; Higgins, M.J.; Farrell, A.; Nakayama, Y.; Jarvis, S.P. Quantitative force 
measurements using frequency modulation atomic force microscopy — Theoretical foundations. 
Nanotechnology 2005, 16, S94-S101. 

3. Adamcik, J.; Berquand, A.; Mezzenga, R. Single-step direct measurement of amyloid fibrils 
stiffness by peak force quantitative nanomechanical atomic force microscopy. Appl. Phys. Lett. 
2011, 98, 193701-193703. 



Sensors 2014, 14 



9378 



4. Xu, X.; Carrasco, C; de Pablo, P.J.; Gomez-Herrero, J.; Raman, A. Unmasking imaging forces on 
soft biological samples in liquids when using dynamic atomic force microscopy: A case study on 
viral capsids. Biophys. J. 2008, 95, 2520-2528. 

5. Melcher, J.; Carrasco, C; Xu, X.; Carrascosa, J.L.; Gomez-Herrero, J.; Jose de Pablo, P.; 
Raman, A. Origins of phase contrast in the atomic force microscope in liquids. PNAS 2009, 106, 
13655-13660. 

6. Kodera, N.; Yamamoto, D.; Ishikawa, R.; Ando, T. Video imaging of walking myosin V by 
high-speed atomic force microscopy. Nature 2010, 468, 72-76. 

7. Shih, P.J. Tip-jump response of an amplitude-modulated atomic force microscope. Sensors 2012, 
12, 6666-6684. 

8. Giessibl, F.J. Forces and frequency shifts in atomic-resolution dynamic-force microscopy. Phys. 
Rev. B 1997, 56, 16010-16015. 

9. Green, CP.; Sader, J.E. Small amplitude oscillations of a thin beam immersed in a viscous fluid 
near a solid surface. Phys. Fluids 2005, 17, doi:l 0.1 063/1. 1995467. 

10. Clough, R.W.; Penzien, J. Dynamics of Structures, 3rd ed.; Computers & Structures, Inc: Berkeley, 
CA, USA, 2003; pp.370-373. 

11. Gotsmann, B.; Anczykowski, B.; Seidel, C; Fuchs, H. Determination of tip-sample interaction 
forces from measured dynamic force spectroscopy curves. Appl. Surf. Sci. 1999, 140, 314-319. 

12. Albrecht, T.R.; Gutter, P.; Home, D.; Rugar, D. Frequency modulation detection using 
high-Q cantilevers for enhanced force microscope sensitivity. J. Appl. Phys. 1991, 69, 

doi: 10. 1063/1. 347347. 

13. Giessibl, F.J. Atomic resolution of the silicon (1 ll)-(7 x 7) surface by atomic force microscopy. 
Science 1995, 267, 68-71. 

14. Hembacher, S.; Giessibl, F.J.; Mannhart, J. Force microscopy with light-atom probes. Science 
2004, 305, 380-383. 

15. Welker, J.; Giessibl, F.J. Revealing the angular symmetry of chemical bonds by atomic force 
microscopy. Science 2012, 336, 444-449. 

16. Giessibl, F.J. A direct method to calculate tip-sample forces from frequency shifts in 
frequency-modulation atomic force microscopy. Appl. Phys. Lett. 2001, 78, doi:10. 1063/1. 1335546. 

17. Tuck, E.O. Calculation of unsteady flows due to small motions of cylinders in a viscous fluid. 
J. Eng. Math. 1969, 3, 29-44. 

18. Sader, J.E. Frequency response of cantilever beams immersed in viscous fluids with applications to 
the atomic force microscope. J. Appl. Phys. 1998, 84, doi: 10. 1063/1. 368002. 

19. Van Eysden, C.A.; Sader, J.E. Resonant frequencies of a rectangular cantilever beam immersed in a 
fluid. J. Appl. Phys. 2006, 100, doi: 10.1063/1.2401053. 

20. Naik, T.; Longmire, E.K.; Mantell, S.C. Dynamic response of a cantilever in liquid near a solid 
wall. Sens. Actuators A 2003, 102, 240-254. 

21. Tung, R.C.; Jana, A.; Raman, A. Hydrodynamic loading of microcantilevers oscillating near rigid 
walls. J. Appl. Phys. 2008, 104, doi: 10.1063/1.3033499. 

22. Riitzel, S.; Lee, S.I.; Raman, A. Nonlinear dynamics of atomic-force-microscope probes driven in 
Lennard- Jones potentials. Proc. R. Soc. Lond. A 2003, 459, 1925-1948. 



Sensors 2014, 14 



9379 



23. Kutana, A.; Giapis, K.P.; Chen, J.Y.; Collier, CP. Amplitude response of single-wall carbon 
nanotube probes during tapping mode atomic force microscopy: Modeling and experiment. Nano 
Lett. 2006, 6, 1669-1673. 

© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.Org/licenses/by/3.0/). 



